Dephasing of a superconducting qubit induced by photon noise 
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We have studied the dephasing of a superconducting flux-qubit coupled to a DC-SQUID based 
oscillator. By varying the bias conditions of both circuits we were able to tune their effective 
coupling strength. This allowed us to measure the effect of such a controllable and well-characterized 
environment on the qubit coherence. We can quantitatively account for our data with a simple 
model in which thermal fluctuations of the photon number in the oscillator are the limiting factor. 
In particular, we observe a strong reduction of the dephasing rate whenever the coupling is tuned 
to zero. At the optimal point we find a large spin-echo decay time of 4fis. 
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Retaining quantum coherence is a central requirement 
in quantum information processing. Solid-state qubits, 
including superconducting ones pi 0, 0, couple to en- 
vironmental degrees of freedom that potentially lead to 
dephasing. This dephasing is commonly associated with 
low- frequency noise 0]. However, resonant modes at 
higher frequencies are harmful as well. In resonance with 
the qubit transition they favor energy relaxation. Off- 
resonance they may cause pure dephasing, due to fluctu- 
ations of the photon number stored in the oscillator. Ex- 
perimentally we show that the quantum coherence of our 
superconducting flux-qubit coupled to a DC-SQUID os- 
cillator is limited by the oscillator thermal photon noise. 
By tuning the qubit and SQUID bias conditions we can 
suppress the influence of photon noise, and we measure 
a strong enhancement of the spin-echo decay time from 
about 100ns to 4/is. 

In our experiment, a flux-qubit of energy splitting hv q 
is coupled to a harmonic oscillator of frequency v v which 
consists of a DC-SQUID and a shunt capacitor [o|.|(J. The 
oscillator is weakly damped with a rate k and is detuned 
from the qubit frequency. In this dispersive regime, the 
presence of n photons in the oscillator induces a qubit 
frequency shift following v q ^ n — v q $ = nSi>a, where the 
shift per photon 8v$ depends on the effective oscillator- 
qubit coupling. Any fluctuation in n thus causes dephas- 
ing. Taking the oscillator to be thermally excited at a 
temperature T and assuming the pure dephasing time 

>> 1/k, we find Q after a reasoning similar to @ : 
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with the average photon number stored in the oscillator 
n = (exphi/p/kT — l) -1 . We note that a similar effect 
was observed in a recent experiment on a charge qubit 
coupled to a slightly detuned waveguide resonator 0. 



When driving the oscillator to perform the readout, the 
authors observed a shift and a broadening of the qubit 
resonance due to the AC-Stark shift and to photon shot 
noise, well-known in atomic cavity quantum electrody- 
namics |lOj . In our experiments, the oscillator is not 
driven but thermally excited. In addition, we are able to 
tune in-situ the coupling constant and Svo, and therefore 
to directly monitor the decohering effect of the circuit. 

Our flux-qubit consists of a micron-size supercon- 
ducting aluminum loop intersected by four Josephson 
junctions 0, 0] fabricated by standard electron-beam 
lithography and shadow evaporation techniques (see fig- 
ure ; note that compared to earlier designs , we 
added a fourth junction to restore the qubit-SQUID cou- 
pling symmetry [lj). When the magnetic flux thread- 
ing the loop Q x sets the total phase across the junc- 
tions 7 g close to 7r, the loop has two low-energy eigen- 
states, ground state |0) and excited state |1) 0, 13- 
The flux-qubit is characterized by the minimum energy 
separation hA between |0) and |1), and the persistent 
current I p ^lj- I n the basis of the energy eigenstates 
at the bias point 7 g = tt, the qubit hamiltonian reads 
H q = -(h/2)(Aa z +ea x ), where e = (I p /e)( lq - n)/{2n) . 
The energy separation is E\ — Eq = hv q = h\J A 2 + e 2 . 
Note that dv q /de = when the qubit is biased at e = 
so that it is to first order insensitive to noise in e, in 
particular to noise in the flux $ x . This is similar to the 
doubly optimal point demonstrated in the quantronium 
experiment |2|. 

The qubit is inductively coupled to a SQUID detec- 
tor with a mutual inductance M (large loop in figure 
and to an on-chip antenna allowing us to apply 
microwave pulses. The readout scheme and the exper- 
imental setup have been described elsewhere 3J. The 
average persistent current in the qubit loop, with a sign 
depending on its state \k) (k = 0,1), generates a flux 
which modifies its critical current Ic ~ lfJ,A to a value 
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FIG. 1: (a) SEM picture of the sample. The flux qubit is the 
small loop containing four Josephson junctions in a row ; the 
SQUID is constituted by the outer loop containing the two 
large junctions. The bar equals 1/im. (b) Measuring circuit 
diagram. The SQUID, represented by its Josephson induc- 
tance Lj, is shunted by an on-chip capacitor C s h through 
superconducting lines of inductance L, forming the plasma 
mode. 



I c ; a bias current pulse of amplitude I m chosen so that 

Iq < I m < Iq allows us to discriminate between the 
two states by detecting the switching of the SQUID. Be- 
fore the measurement, when the bias current It, < Ic, 
the SQUID behaves as a Josephson inductance Lj which 
depends on the flux threading it and on lb- It is con- 
nected to an on-chip capacitor C s h through a line with a 
stray inductance L (see figure ^d) and thus forms a har- 
monic oscillator of frequency v p — l/2wy/ (L + Lj)C s h 
called the plasma mode 0, |(j (note that the junction 
capacitance is much smaller than C s h)- We can write 
its hamiltonian H p = hv p a^a, where a (a^) is the anni- 
hilation (creation) operator. The total current flowing 
through the SQUID is thus h + i, with i = 5io(a + a') 
being the operator for the current in the plasma mode 
and Sio the rms fluctuations of the current in the oscil- 
lator ground state 5io — \Jhv p j2{L + Lj). The SQUID 
circuit is connected to the output voltage of our wave- 
form generator E via an impedance Z in , and to the in- 
put of a room-temperature amplifier through Z out which 
define the oscillator quality factor Q = 2ttv p /k. Zi n 
and Zn Ut take into account low-temperature low-pass fil- 
ters 0], and on-chip thin-film gold resistors ther- 
malized by massive heat-sinks. The resulting impedance 
seen from the plasma mode is estimated to be 9ktt at 
low frequencies and of order 50017 at GHz frequencies. 
The measurements were performed at a base temperature 
T b = 30mK. 

The applied magnetic field and the bias current A re- 
sult in a circulating current J in the SQUID loop [Tj|- 
Via the qubit-SQUID coupling M the qubit phase j q 
will be affected, so that we can write the qubit en- 
ergy bias as a sum of two contributions e = r\ + A, 
where r\ = 2I p (Q x — $o/2)//i is controlled by Q x and 
A = 2I p MJ(Ib)/h only depends on if, [l^. This depen- 
dence has two important consequences. First, the qubit 
bias point is shifted by the measurement pulse, allow- 
ing us to operate the qubit at the flux-noise insensitive 
point while keeping a measurable signal Second, it 
gives rise to a coupling between the qubit and the plasma 
mode described by a hamiltonian Hj = h[gi(Ib)(a+<r) + 
g2(h){a + a^) 2 ]a x , where gi{h) = (l/2)(dX/dI b )Si , and 




FIG. 2: (a) Top : Principle of the spectroscopy experiments : 
a bias current pulse of amplitude Ib p i < Ic ~ \\iA is applied 
while a microwave pulse (MW) probes the qubit resonance 
frequency. The qubit state is finally measured by a short 
bias current pulse as discussed in Q. Bottom : Qubit spec- 
troscopy for I^i varying between 0/j,A to 0.6/j.A with steps of 
0.1/j.A (bottom to top). The curves were offset by 100MHz 
for clarity. The solid curves are fits to the formula for v q . 
(b) Curve X(h) deduced from the spectroscopy curves as ex- 
plained in the text. The solid line is a parabolic fit to the data. 
The decoupling condition is satisfied at I£ = 180 ± 20nA. (c) 
Calculated frequency shift Svo(Ib, e) for the parameters of our 
sample. The white scale corresponds to —20MHz, the black 
to +40MHz. Along the dotted line e m (h), Su = 0. 

.92(4) = {l/A){d 2 X/dll){5i ) 2 7\. We note that this cou- 
pling hamiltonian depends on 1^ via <?i and gi and is thus 
tunable in-situ. In particular it is possible to cancel g\ by 
biasing the SQUID at a current 1% such that dX/dlb = 0. 
The qubit is then effectively decoupled from its measuring 
circuit 13]. Our design therefore allows us to study the 
effect of the coupling between the qubit and its measuring 
circuit by varying 7^, while keeping all other parameters 
unchanged. 

To obtain the coupling constants g\ and g2, we per- 
formed extensive spectroscopic measurements of the 
qubit, as a function of both If, and $ x . We applied 
a pre-bias current pulse Ib P i through the SQUID while 
sending a long microwave pulse, followed by a regular 
measurement pulse : 3J at a value I m (see figure^). We 
measured the SQUID switching probability as a func- 
tion of the microwave frequency, and recorded the po- 
sition of the qubit resonance as a function of Ib P i and 
Q x . The data are shown in figure [2t for various val- 
ues of Ib p i- We observe that for each bias current, a 
specific value of external flux & x (Ib p i) realizes the op- 
timal point condition. Fitting all the curves with the 
formula v q = ^ A 2 + {MJbpT) + 2I P {$ X - $ /2)//i] 2 , we 
obtain the qubit parameters M = 6.5pH, A = 5.5GHz, 
I p = 240nA, and also A(Jfc) which is shown in figure 
13) together with a parabolic fit. Decoupling occurs at 
lb = 180 ± 20nA and not at lb = because of a 4% 
asymmetry of the SQUID junctions. We also measured 
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FIG. 3: (a) Qubit line shape at the optimal point. The solid 
line is a fit assuming a double lorentzian shape, (b) Rabi 
oscillations (frequency 100MHz) at the optimal point. The 
inset shows well-behaved oscillations with nearly no damping 
during the first 100ns. (c) Measurement of Ti at the optimal 
point ; the dashed grey line is an exponential fit of a time 
constant 4fis. (d) Spin-echo pulse sequence and signal at the 
optimal point. 



the parameters of the SQUID oscillator by performing 
resonant activation measurements and fitting the depen- 
dence of the resonant activation peak as a function of 
h and § x @. We found a maximum plasma frequency 
v p = 3.17GHz, C sh = 7.5 ± 2pF and L = 100 ± 20pH, 
consistent with design values. The width of the peak 
also gives us an estimate for the oscillator quality factor, 
Q = 120±30. 

From the previous measurements we know the param- 
eters of the total hamiltonian H — H q + H p + Hj and 
we can deduce the value of Svq by second-order per- 
turbation theory : 8v$ = 4[(gi(Jb) sinfl) 2 ^ a — 

<72(^h)cos#] where 9 is the mixing angle, defined by 



cos( 



e/Ve 2 + A 2 . The first term in 5vq is the usual 



AC-Zeeman shift obtained without using the rotating 
wave approximation which is not valid in our case. Note 
in particular that the sign of this shift only depends on 
the sign of v q — v v which in our experiments is always pos- 
itive. The second term is due to the dependence of the 
SQUID Joscphson inductance on the qubit state [jj 0] , 
and it has the same sign as e since gi is negative. There- 
fore, for some value e m (Ib) < 0, one obtains Svq = 0. 
This is shown as a dashed line in figure [2t in which we 
plot 5vo(e,Ib)- If dephasing is indeed limited by ther- 
mal fluctuations, we expect the dephasing time to be 
maximal along e m (I&). We note that the curve includes 
(7(, = /^e = 0), so that this bias point is optimal with 
regard to bias current noise, flux noise and photon noise. 

We now turn to the measurements of the qubit coher- 
ence properties around this optimal point, as character- 
ized by the relaxation time T\, the qubit spectral line 
shape, and the spin-echo decay time T ec h |l9j |. The line 
shape was measured using a long microwave pulse (2[is) 



at a power well below saturation. Figure^ shows a typ- 
ical result at the optimal point. For this specific sample, 
we observed a twin peak structure which likely results 
from one strongly coupled microscopic fluctuator. In ad- 
dition, the width of the line as well as the average value of 
the gap A changed significantly in time, which indicates 
that the residual linewidth is probably due to a larger 
number of fluctuators more weakly coupled. We stress 
that we observed the splitting all along the v q (<& x ) spec- 
trum in contrast to [lq. Fitting the peaks to a sum of 
two Lorentzians of widths w\ and wi we define an effec- 
tive dephasing time t% = 2/n(wi + w%). At the optimal 
point t% varied between 50 and 200ns. 

Despite the fluctuators, we were able to induce Rabi 
oscillations by applying microwave pulses at the mid- 
dle frequency of the split line. An example is shown 
in figure |3Jd at the optimal point. The oscillations decay 
non-exponentially and display a clear beating. Neverthe- 
less, by driving the qubit strongly enough, we could ob- 
serve well-behaved oscillations for hundreds of nanosec- 
onds (see inset of figure ■ We measured the energy 
relaxation time T\ by applying a 7r pulse followed after a 
delay Dt by a measurement pulse (see figure 0). At the 
optimal point, we found that T% — Afis. To quantify the 
dephasing further we also applied the spin-echo sequence 
[l9|. depicted in figure^. Spin-echo measurements are 
particularly relevant for our purpose, because the pho- 
ton noise in the plasma mode occurs at a relatively high 
frequency set by k ~ 130MHz. In such conditions, this 
noise affects the spin-echo damping time T ec ho as strongly 
as Ramsey experiments so that T ec ho is also given by for- 
mula n 7] ; on the other hand the spin-echo experiment 
is not sensitive to the low-frequency noise responsible for 
the qubit line splitting. The results are shown in fig- 
ure 01 at the optimal point, by a set of curves obtained 
at different delays At between the it pulse and the last 
7r/2 pulses. Fitting the decay of the echo amplitude as a 
function of the delay between the two 7r/2 pulses with an 
exponential, we find T ec ) lo — 3.9±0.1/xs. Compared with 
previous experiments on flux-qubits the long Rabi 
and spin-echo times were obtained by reducing the mu- 
tual inductance M, and biasing the qubit at the optimal 
point. 

We studied the variation of T\ as a function of the bias 
current at the flux-insensitive point e = 0. This required 
us to adjust the flux at the value & x (Ib). Results are 
shown in figure^,. We observed a clear maximum of Ti 
for Ii, = lb - This demonstrates that at least part of the 
qubit relaxation occurs by dissipation in the measuring 
circuit. We then investigated the dependence of T ec ho 
and t% on e for lb = Ij* (figure^ top, full circles and full 
squares). As expected, we observe a sharp maximum for 
T ec ho at e = and a shallow one for ti. However, at a 
different bias current lb — fj,A, the maximum of T ec ho 
and ti is clearly shifted towards e < 0. We measured 
the position of this maximum in ti as a function of lb as 
shown in solid squares in figure For dephasing caused 
by flux noise or bias current noise, the maximal coherence 
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FIG. 4: (a) Measurement of Ti versus h at the flux-noise 
insensitive point e = 0. (b) Measurement of T ec ho (circles), ti 
(squares) and of the qubit frequency (triangles), as a function 
of e for It = II (top) and lb = 0/j.A (bottom). The dotted 
line is a fit to the formula for v q ; the solid black line is the 
prediction of equation for T = 70mK and Q = 150. (c) 
Value of e for which ti is maximum (full squares) compared 
to the theoretical e m (I b ) (full line). 

time should always be obtained at e = ; the observed 
deviation proves that a different noise source is active in 



our experiments. We find that thermally induced photon 
number fluctuations in the plasma mode explains our re- 
sults. In figure 0J: we draw the curve e m (h) of figure 
|2Ji, where the photon-induced shift 5vq equals (solid 
line). The agreement between the data points and this 
curve, obtained directly from measured parameters, is 
excellent. In addition, assuming a reasonable effective 
oscillator temperature of T — 70mK |5j and a quality 
factor of Q = 150, which yields a mean photon number 
n = 0.15, the dephasing time predicted by equation 
n closely matches the spin-echo measurements both for 
h = II and h — [iA [l7| (see the solid line in figure^). 
We stress that even at such small n the photon number 
fluctuations can strongly limit the qubit coherence. This 
suggests that increasing the plasma frequency could lead 
to significant improvement. 

In conclusion, we present experimental evidence that 
the dephasing times measured in a flux-qubit can be lim- 
ited by thermal fluctuations of the photon number in the 
SQUID detector plasma mode to which it is strongly cou- 
pled. By careful tuning of flux and current bias, we could 
decouple the qubit from its detector and reach long relax- 
ation and spin-echo damping times (4/is). These results 
indicate that long coherence times can be achieved with 
flux qubits. 
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